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Abstract 



A conserved cosmological perturbation is associated with each quantity whose 
local evolution is determined entirely by the local expansion of the Universe. It 
may be defined as the appropriately normalised perturbation of the quantity 
defined using a slicing of spacetime such that the expansion between slices is 
spatially homogeneous. To first order, on super-horizon scales, the slicing with 
unperturbed intrinsic curvature has this property. A general construction is 
given for conserved quantities, yielding the curvature perturbation £ as well 
as more recently-considered conserved perturbations. The construction may 
be extended to higher orders in perturbation theory and even into the non- 
perturbative regime. 
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I. INTRODUCTION 



Observation of the peak structure in the Cosmic Microwave Background (CMB) 
anisotropy has now confirmed that cosmological perturbations are present before the rel- 
evant scales enter the horizon, with an almost flat (scale-invariant) spectrum [1-3]. The 
only known explanation for this state of affairs is that the perturbations originate during an 
almost exponential inflation, from the vacuum fluctuation of one or more light scalar fields. 1 
In the simplest case only one light field is responsible for the perturbations observed, either 
the inflaton or some other field . 

According to this explanation, classical cosmological perturbations first come into exis- 
tence a few Hubble times after horizon exit during inflation. At that stage the situation 
is very simple; each light field (defined as one with an effective mass much less than the 
Hubble parameter H) has a Gaussian perturbation with an almost flat spectrum, (H/2n) 2 . 
The problem is to evolve this simple initial condition forward in time to the primordial nu- 
cleosynthesis epoch, in the face of our ignorance about the detailed evolution of the Universe 
before nucleosynthesis. 

Fortunately, scales of cosmological interest are still far outside the horizon at nucleosyn- 
thesis. As a result there exist perturbations which are under suitable conditions conserved, 
and largely avoiding the need for more detailed information. One of these [11-14] is the 
'curvature perturbation' (, which is associated with the perturbation in the total energy 
density p. 2 In the usual case that ( originates from the perturbation in the inflaton field, it 
is supposed to be conserved between the end of inflation and the primordial era, and in the 
alternative curvaton scenario [16,17] (see also [18,19]) ( is supposed to be conserved after 
the curvaton decays. 3 Recently, further conserved quantities Q and Q have been considered, 
that are associated with the perturbations in individual energy densities pi [14] and number 
densities rij of conserved quantities [23]. The conservation of the former is invoked in the 
curvaton scenario, during the era when the curvaton field is oscillating and ( is growing. 
The latter are invoked when considering possible isocurvature components of the primordial 
density perturbation. 

In this paper, we present a unified treatment of the conserved quantities (, Q and Q, which 



X A related hypothesis replaces inflation by an era of collapse ('pre-big-bang' [4-6] or 'ekpyrotic' 
[7-9]), but there is so far no accepted theory of a bounce and therefore no firm prediction from 
collapsing cosmologies. In particular, there is so far no accepted string-theoretic description of a 
bounce [10]. 

2 The quantity C defines the curvature perturbation on spacetime slices of uniform energy density. 
As we discuss in Section V, on super-horizon scales it is practically the same as 1Z which defines 
the curvature perturbation on slices orthogonal to comoving worldlines. The latter quantity is the 

0m Of [15]. 

3 An analogous scenario has been proposed in the pre big bang scenario [20,21]. In this scenario 
though, the required scale-invariant curvaton field perturbations will be generated only if the 
curvaton has a non-trivial coupling and for particular initial conditions [22,6]. 
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is more complete than anything that has been given before. Taking the particular example 
of ( as a starting point, we begin in Section II by showing how, to any order in cosmological 
perturbation theory, conserved quantities may be constructed from perturbations that are 
defined on a spacetime slicing of uniform integrated expansion. In Section III we generalise 
the construction and consider the conserved quantities Q and Q. In Section IV, we show 
that in the usual case of first-order perturbation theory, the spatially flat slicing is one of 
uniform expansion if the shear of the worldlines is negligible. In Section V we consider the 
comoving shear, and show that it is expected to be negligible in the entire super-horizon 
regime. We conclude in Section VI. An Appendix mentions some peripheral issues. 



II. ENERGY CONSERVATION AND THE CURVATURE PERTURBATION 

In this section we explain the general principle which allows us to construct conserved 
quantities. We focus on the particularly important example of the curvature perturbation 
( [11-14], after which it is clear how other conserved perturbations may be constructed. 

The curvature perturbation ( is so called because it defines the curvature perturbation 
on slices of uniform energy density [14]. Equivalently though, via the gauge transformations 
of Section IV, it defines the energy density perturbation on spatially flat slices, according to 
the formula [12] 

(= 5p . (i) 

^ 3(p + P) {) 

This definition is the one that we shall use. 

Our starting point is the energy continuity equation. In an unperturbed Friedmann- 
Robert son- Walker (FRW) universe the continuity equation for the energy density p takes 
the form 

p=-3H(p+P), (2) 

where H is the Hubble expansion rate and P is the pressure. In the real perturbed Universe, 
the same equation (2) still holds along each comoving worldline, so long as the dot is taken 
to denote the derivative with respect to the proper time, r, along the comoving worldline 
and we define H locally through the equation 

where V is an infinitesimal comoving volume. Equivalently, the local continuity equation 
may be written as 

V^ = -(p + P),. (4) 



or 



^ = -3(P+^). (5) 
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where N is the local logarithmic integrated expansion (the number of Hubble times) defined 



as 

N 



J Hdr . (6) 



Our crucial assumption now is that the pressure perturbation is practically adiabatic. 
This assumption means that the local pressure, P, is a practically unique function of local 
energy density, p. i.e., 

p = Hp) , (7) 

where P is the same function for all worldlines. This allows Eq. (5) to be integrated. Setting 
N = on an initial spacetime slice the integration gives p as a unique function of the local 
integrated expansion, N, upto an initial integration constant 

p = p(N + 5N), (8) 

where the integration constant for each worldline, SN, is determined by the actual density 
on the initial hypersurface, p|at=o = p(SN). 

Subsequent spacetime slices of fixed N correspond to a uniform integrated expansion slic- 
ing of the spacetime 4 , meaning that the integrated expansion going from one slice to another 
is spatially homogeneous. For linear perturbations about an FRW cosmology, there is an 
infinity of such uniform- N slicings, since we can start with any initial slice and propagate it 
by calculating N from that slice along each comoving worldline. In Sections IV and V, we 
show that on super-horizon scales a particular uniform- N slicing is the uniform curvature 
slicing (i.e., the one with unperturbed intrinsic scalar curvature). In what follows we will 
restrict our attention to spatially flat FRW models and will refer to this as the spatially flat 
slicing. 

Now we come to the crucial point. When evaluating the density, p, on any uniform- N 
slicing, the perturbation 5N of the quantity appearing in Eq. (8) is time-independent, by 
construction. This statement holds to any order in cosmological perturbation theory so long 
as one can construct a uniform- iV slicing along the comoving worldlines. 

Writing 5N in terms of the density perturbation on spatially flat slices, to first order, 
one finds the conserved quantity 

dp r p'(N) p 

which is — C defined in Eq. (1). To arrive at the conserved quantity (, we considered the 
flat slicing. Were we instead to use some other uniform-expansion slicing, the conserved 
quantity defined by the right hand side of Eq. (9) would be different from (, but it would be 
related to ( by the gauge transformation Eq. (40). Hence it would be conserved if and only 



4 Note that this is not the same as the uniform Hubble slicing introduced by Bardeen [15,11] which 
refers to the local expansion rate of the normals. 
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if ( is conserved, and we lose no generality by fixing the choice of the uniform-expansion 
slicing as the flat one. 

The constancy of ( (on sufficiently large scales and assuming that the pressure pertur- 
bation is adiabatic) was obtained several years ago [11] in the context of Einstein gravity. 
More recently, its constancy under the same condition was obtained directly from the local 
conservation of energy [14] using a purely geometric argument equivalent to the one that we 
have given. In the present paper, we are going to show in Sections IV and V that in this 
context all super-horizon scales are 'sufficiently large'; in other words, we will show that the 
flat slicing is a uniform-expansion one on all super-horizon scales. 

It is worth noting that the conservation of ( can hold in even more general circumstances, 
because it comes from the generalised adiabatic condition Eq. (8) which may hold even if 
the energy conservation equation (5) fails. Thus, ( will be conserved even if there is an 
additional source term Q on the right-hand-side of in Eq. (energy), so long as Q (the energy 
transfer per Hubble time) is itself a unique function of the local density for all worldlines, 
i.e., Q = Q(p), as reported in Ref. [24]. 

Going to second order, and again working on some uniform-iV slicing, the conserved 
quantity is 



This second-order extension of the conserved quantity ( has not been given before. It will be 
useful in propagating forward the evolution of second-order perturbations produced during 
inflation [25-27] through the end of inflation and relating them to observations. Also, we 
note that Sasaki and Tanaka [25] have shown that it is possible to use a uniform- A slicing 
to study non-linear field perturbations on large scales during inflation. 



III. OTHER CONSERVED QUANTITIES 

Generalising from the construction of (, it is clear that for any monotonically increasing 
or decreasing quantity, satisfying a local conservation equation of the form 

V§£ = !/(/)> (12) 
we can construct a conserved first-order perturbation 

X f = -H 5 -l, (13) 

with 5f evaluated on some uniform- A slicing which we will take to be the spatially flat one. 
This construction gives ( = X p as a special case, and we shall now see how it gives the other 
conserved quantities Q and Q [14,23]. 
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A. Separately conserved energy densities 



Suppose the total energy density p of the Universe is a sum of components, p i: each one 
of them either radiation or matter, and with no energy transfer between the components. 
In that case the pressure of each component is a unique function of its energy (Pj = p^/3 
for radiation and p = for matter) and each component satisfies its own separate energy 
conservation equation (since there is no energy transfer) 5 

V^ = -(P, + P). (14) 
As a result there are the separately conserved perturbations 



Q = X n (15) 

Pi 

5p 



-H*» (16) 

Pi 



3pi + Pi ' 



(17) 



where 5pi is evaluated on the flat slicing. This is another result of [14]. 

One can express the total density perturbation, (, as a weighted sum of the separate Q; 

C = ^- (18) 

l^i Pi 

If the Q are all equal, then C — C« which is constant. Otherwise £ may have some variation, 
determined by the conserved isocurvature perturbations defined by 

^. = 3(0-0)- (19) 

The condition that the Q are equal is just the adiabatic condition, that all of the separate 
energy densities (and hence the total pressure) are uniform on slices of uniform total energy 
density. 

There are two eras in the early Universe where separately-conserved Q have been invoked. 
One is the comparatively late era, beginning when the temperature falls below 1 MeV and 
ending when cosmological scales start to approach the horizon. 6 

The energy density during this era has four components, 

P = Pcdm + Pb + Pv + P 7 , (20) 



5 In this expression, V is the volume which is comoving with the flow of pi. This is not strictly 
the same as the volume which is comoving with the flow of total energy density, but we are going 
to show in Section IV that on super-horizon scales all comoving volumes become equivalent. 

6 Recall that electron-positron annihilation and neutrino decoupling both take place when the 
temperature is around 1 MeV. 
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with the radiation (photons and neutrinos) dominating the matter (Cold Dark Matter and 
baryonic matter). The values of the four conserved quantities Ccdm, Cb, Cu and £ 7 determine 
the evolution of the entire set of cosmological perturbations after horizon entry, and can 
therefore be determined by observation. The three isocurvature perturbations (convention- 
ally defined relative to the photon density to be S c d m = 5 c dm7, $b = Sb-/ and S v = S U1 ) 
are found by observation to be at most of order ( [29-31,2]. Since radiation dominates, one 
deduces from Eq. (18) that ( is constant on large scales to high accuracy during this era. 

The other era, which occurs in the recently proposed curvaton scenario [16,17] (see also 
[18,19]), is the era (after inflation, but before primordial nucleosynthesis) when the massive 
curvaton field, a, oscillates (P CT = 0) in a radiation background (P r = p r /3), 

P = Pa + Pr • (21) 

Here, ( T is supposed to be negligible so that the total curvature perturbation, (, is given by 
Eq. (18) as 

= z^zV c - • (22) 

Well before the curvaton decays, the radiation is supposed to dominate so that ( grows like 
Pa/Pr oc a(t), providing an example where the total ( is not conserved on super-horizon 
scales after inflation. 



B. Conserved number densities 

If rii is a conserved number density, then rij is inversely proportional to the volume. 7 
A conservation law of the form given in Eq. (53) is satisfied with y(f) = —f and / = n«, 
leading to the conservation of the first-order perturbation [23] 

C.^„, = i^. (23) 

These conserved quantities find an application [23,32] in connection with the three isocur- 
vature perturbations iS c d m) <~>b and S v . Before T ~ 1 MeV, these quantities are not the 
appropriate ones to consider, because the separate energy density perturbations Q may vary 
with time or be simply undefined. (The latter is the case for (b before the quark-hadron 
transition.) One can however consider instead the number density n c d m of Cold Dark Matter 
particles, the density % of Baryon Number and the density rii of Lepton Number. Each 
of these number densities corresponds to a conserved quantity after some epoch, which may 
be regarded as the epoch when the quantity originates. The corresponding perturbations Q 



7 The volume should be the one comoving with the flow of the conserved quantity, but as already 
stated we are going to show in Sections IV and V that the choice of comoving volume is irrelevant on 
super-horizon scales. This irrelevance is assumed implicitly when the conservation of Q is discussed 
in [23]. 
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are thus conserved, and after the temperature falls below 1 MeV they determine the three 
isocurvature perturbations according to the formulas [23] 




(24) 



±s B = Cb-C 



(25) 




(26) 



In the last formula, £ is the lepton asymmetry which must satisfy the nucleosynthesis con- 
straint |£| < 0.07. In this way, the three isocurvature perturbations can be calculated (or 
shown to vanish) given a model of the early Universe. 



The main goal of this section is to show that the local expansion of the Universe between 
spatially flat slices is uniform on sufficiently large scales where shear is negligible. This result 
is purely geometric, making no reference to the theory of gravity. Our treatment amplifies 
the original one in Ref. [14], and in particular we show for the first time that the result is 
independent of the spacetime threading that defines the expansion. 



The unperturbed FRW metric and comoving coordinates are defined by the line element 



corresponding to metric components goo — — 1, goi — and gij = 5ijd 2 . 

We are interested in the perturbed spacetime that is our Universe, which we assume can 
be described by linear perturbations about a FRW geometry. To define the perturbations 
one has to choose a coordinate system which reduces to Eq. (27) in the limit where the 
perturbations vanish. Such a coordinate system (gauge) defines a time-slicing (the spatial 
hypersurfaces with constant time coordinate) and a threading (the worldlines with constant 
space coordinates) of the spacetime. Since the coordinate system is required to coincide 
with Eq. (27) in the limit where the perturbations vanish, the slicing and threading coincide 
with the unperturbed ones in that limit. We shall take this requirement for granted when 
referring to a 'generic' slicing or threading. 

Once the perturbations are defined, their evolution to first order may be described using 
the unperturbed coordinate system, and Fourier components with different wavevectors k 
decouple. The super-horizon regime is the regime aH/k 3> 1. 

In this paper we are interested in the scalar mode of the perturbations in the metric (no 
gravitational waves or vorticity). In a generic gauge the Fourier components of the metric 



IV. UNIFORM EXPANSION BETWEEN FLAT SLICES 



A. The metric perturbation 



ds 2 = -dt 2 + a 2 (t)6 ij dx i dx i , 



(27) 
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perturbation are specified by functions A, B, D and E, 

1 



k 

where Pij projects out the traceless part; 



2 5g 00 = -A (28) 

2a- 2 5g 0t = -B t = %k { B (29) 

l -a~ 2 bg i3 = 5ijD + P tj E (30) 

= ^-^-^E, (31) 



P^-^ + ^j, (32) 



and 



-^ = D + ^E. (33) 

Under the coordinate transformation t — > t + At and x l — > x % + Ax*, with the Fourier 
component of Ax 1 of the form 

Ax l = -i^Ax, (34) 
the metric components transform according to 

AA = -At (35) 
AB = aAx + kAt (36) 

AD = -^Ax - HAt (37) 

AE = kAx (38) 
AiP = HAt (39) 

The perturbation E can be eliminated by a transformation of the spatial coordinates, while 
ip depends only on the slicing. The latter defines the intrinsic scalar 3-curvature of the 
slicing, and is called its curvature perturbation. 

The perturbation in a generic quantity g defined in the unperturbed Universe, such as 
energy density, pressure or the value of a scalar field, has the transformation 

Ag = -gAt. (40) 

Applied to the energy-density, this transformation along with Eq. (39) leads to the gauge- 
invariant definition of the 'curvature perturbation' (, 



8 These are the quantities defined in [33], and are equal respectively to the quantities A, B(°\ Hl 
and H^p of Bardeen [15]. The quantities Ax and At below are respectively the the L and aT of 
Bardeen, and the quantity V in Eq. (60) is the of Bardeen. 
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C = -^-#-. (41) 
P 

Evaluated on slices of uniform density it is indeed the curvature perturbation, but evaluated 
on flat slices it specifies instead the energy density perturbation through Eq. (1). 



B. Shear and the expansion rate 

A given threading of spacetime is associated with an expansion 9, a traceless symmetric 
shear and an antisymmetric vorticity ou^. At a given spacetime point, in a locally- inert ial 
rest-frame, these quantities are given by the standard decomposition [33] 

diWj = ^95ij + Oij + Uij, . (42) 

where Wi is the three- velocity of the infinitesimally nearby threads. The expansion 9 gives 
the rate of increase with respect to proper time r of an infinitesimal volume V expanding 
with the threads, 

= = (43) 

For the scalar perturbations that we are considering, the vorticity vanishes, and the shear 
takes the form 

(Tij = P l3 o . (44) 

We shall call a the shear as well. 

In a given gauge, the metric perturbations determine the shear and expansion rate of 
the coordinate threads according to the expressions [34] 

a = E (45) 
86 = 3D- 3HA . (46) 

The second relation may be written as 

59 = -3ip -a- 3HA . (47) 

Going to a new threading corresponding to the transformation Eq. (34), the change in 
the Fourier component of the local velocity field of the threads is 

k- 

Awi = Ax i = Ax , (48) 

The corresponding changes in the shear and the expansion are equal and opposite, 

A(59) = —Act = kAx , (49) 

while the perturbations ip and A are unchanged. It follows that Eq. (47) is valid for any 
threading, not only for the coordinate threading. (In [14], Eq. (47) was given for the special 
case of the threading normal to the coordinate slicing.) 
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Following [14] we consider, instead of 9, the expansion 9 with respect to coordinate time. 
Its perturbation is 



89 = -3ip - a . 



(50) 



For the flat slicing (ip = 0) this becomes 



59 = -a. 



(51) 



This result is true for any choice of the threading that defines the expansion rate, and has 
been derived without any reference to a theory of gravity. 

The expressions given so far are valid quite generally. We are interested though in super- 
horizon scales. On sufficiently large scales the shear must become negligible compared with 
the Hubble parameter, so that we recover an unperturbed FRW universe. It follows that 
89 is negligible on sufficiently large scales, or in other words that the expansion between 
successive flat slices becomes unperturbed. This means that on sufficiently large scales we 
can make the approximation 



where H(t) is the usual unperturbed quantity and 9 is the expansion with respect to co- 
ordinate time of a generic threading. (Remember that we consider only threadings which 
coincide with the unperturbed one in the limit of zero perturbation.) 

Using Eq. (52), we can combine the results of section II and III to derive the following 
general result for cosmological perturbations: 

Consider a monotonically increasing or decreasing quantity /, defined in some 
region of spacetime, and its first-order perturbation Sf defined on the spatially 
flat slicing. Consider also some threading of spacetime, defining an infinitesimal 
volume element V. If / satisfies a local conservation equation of the form 



0(x,t) = 3H(t), 



(52) 




(53) 



then the rate of change of the perturbation 




(54) 



is 




1 



(55) 



where a is the shear of the threading. As a result, Xf is conserved on sufficiently 
large scales, where the shear is negligible. 
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The variation of ( 



Taking the derivative of Eq. (41) and using the local conservation of energy along co- 
moving worldlines one finds 

C = -^^P nad -a, (56) 

where a is the shear of the comoving worldlines and the non-adiabatic part of the pressure 
perturbation is 

SP nad = 5P--5p. (57) 
P 

This result was derived by essentially the above method in [14]. It was first derived (by 
a different method, and actually for the curvature perturbation 1Z) in [15], Eqs. (5.19) to 
(5.21). In the particular case that the Universe consists entirely of matter and radiation, ( 
is given by Eq. (18). One easily checks [23] that this expression is compatible with Eq. (56). 

V. SHEAR ON SUPER-HORIZON SCALES 

According to Eq. (55), the quantity Xf is conserved on scales which are sufficiently large 
that a is negligible. In this section we argue that any super-horizon scale is 'sufficiently 
large' in this context. To be more precise, we argue that the shear satisfies 

\a\/H <^{k/aH) , (58) 

which ensures that in one Hubble time the change in Xf is less than k/aH. 9 

In making the argument, we shall invoke the Einstein field equations. This is not much 
of a restriction as we do not necessarily have to specify any physical origin for the stress- 
energy tensor, but simply equate it with a fixed multiple of the Einstein tensor derived 
from the metric. This leads to a purely geometrical definition of, e.g., the comoving density 
perturbation, which need have nothing to do with the motion of particles. However at 
relatively late cosmic times, it may be safe to assume that the Einstein field equations are 
satisfied with the stress-energy tensor related to particle physics content in the usual way. 

The other assumption we make is that anisotropic stress is negligible on super-horizon 
scales. As noted by Bardeen in 1980 [15], significant anisotropic stress on super-horizon scales 
would generate shear and the curvature perturbation, but there is no known mechanism for 
generating such stress. 

We first note that any spatial gauge transformation, Eq. (34) corresponds to a change 
in the local physical velocity Avi = aAwi (Eq. (48)). This generates a change in the 



9 On the left hand side of this expression, a is the typical magnitude of the shear on scale k, 
defined for instance as 10 the square root of its spectrum V a {k). 
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shear \Aa\/H = vk/aH). Since we are dealing with small perturbations, C 1 so that 
\Aa\/H < k/aH. 

It follows that we need only establish Eq. (58) for the shear of the comoving threading, 
which a shall denote from now on. Generalising the discussion of Bardeen [15] to include 
the case where Pj p may vary, we shall show that in fact 

(k/aH) 2 . (59) 

The comoving shear is related to two commonly used gauge-invariant variables, namely 
the curvature perturbation —TZ of slices orthogonal to comoving worldlines (comoving slices) 
and the curvature perturbation $ of zero-shear hypersurfaces (the Bardeen potential): 11 

(The quantity V defines the velocity Vi of the comoving worldlines relative to the the zero- 
shear threading, through the relation Vi = —i{ki/k)V [33].) The curvature perturbation TZ 
is closely related to the curvature of uniform- density slices, (: 

p 

where the subscript 'com' denotes the comoving slicing. 

The combined energy and momentum constraints of Einstein's equations relate the co- 
moving density perturbation to the Bardeen potential 

HS P™ = 2 (±)\, (62) 
p 9(1 + w) \aH) ' K ' 

where w = P/p. Thus we can rewrite Eq. (60) for the comoving shear in terms of ( and the 
Bardeen potential, giving 



a 
H 




k V 



9(1 + w) \aH / 



$ . (63) 



Equation (63) ensures that the comoving shear will be small (<j/H <C () on super-horizon 
scales so long as the Bardeen potential, <3>, remains of the same order as (. Notice also that 
if the <3> remains finite on large scales, then the comoving density perturbation Eq. (62) also 
vanishes on large scales, and TZ and ( are equal (and constant) on large scales. However, it 
has been argued [15] that cosmological perturbation theory can still be valid even if certain 
curvature perturbations, in particular <3>, become formally bigger than one. 

The Bardeen potential is not uniquely determined by the value of (, but the Einstein's 
equations give a first-order evolution equation [33] (in the absence of anisotropic stress [15]) 



n We are defining —TZ = ip with the right hand side evaluated on comoving slices, which corre- 
sponds to established conventions. 
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H- 1 $ + 



5 + 3w 




$ = --(i + «,)C. 



(64) 



During conventional slow-roll inflation (with w ~ — 1 and $ <C the Bardeen potential 

is indeed small, $ ~ —3(1 + w)(/2 on super-horizon scales. But we wish to eliminate the 
possibility that the Bardeen potential subsequently becomes large on super-horizon scales. 



1. Adiabatic perturbations 

For strictly adiabatic matter perturbations, with <5P na d = 0, Eqs. (56), (63) and (64) 
yield coupled first-order equations for the evolution of ( and $, which to lowest order in 
k/aH gives 

(65) 

+ = -§(!+»)<, (66) 

which yield two independent long-wavelength solutions, which are represented by 

C^C + , (67) 
$ ~ C -e- {b+ ^ N l 2 , (68) 

where wN = J wdN. The first of these solutions, with constant ( on large scales, remains 
the "growing mode" solution so long as 

H-'C oc c_Jfe 2 e 3 ( 2w -*- 1 W 2 , (69) 

for the "decaying mode" on large scales, approaches zero. This is always true in an expanding 
universe (N — > +oo) so long as w — > < 1, i.e., P < p. This is easily interpreted as 
the condition for the decay of the shear relative to the Hubble rate (cr/H) in an expanding 
universe. 

Using these same equations, we can understand the super-horizon evolution of the shear 
and the curvature perturbations in a collapsing universe (N — > — oo). For w < 1, the 
shear grows relative to the Hubble rate, and ( does not remain constant. The critical case 
w — 1 (maximally stiff fluid) occurs if the energy density is dominated by scalar fields with 
negligible potential. This is supposed to happen in the pre big bang scenario, and in the late 
stages of the second version of the ekpyrotic scenario [8] where the bounce is supposed to be 
singular from the four-dimensional viewpoint. For w — 1, £ on super-horizon scales grows 
logarithmically with respect to cosmic time and has a strongly scale-dependent spectrum 
[5]. (It is however [5] still small at the string epoch, which in the pre big bang scenario is 
supposed to be the bounce epoch.) 

In the first version of the ekpyrotic scenario [7], where the bounce is supposed to be 
non-singular from the four-dimensional viewpoint, collapse is driven by a scalar field with 
a steep negative potential which violates the dominant energy condition and gives w ^> 1. 
The same is supposed to happen in the second version of the ekpyrotic scenario [8] at early 
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times. In these cases, the shear rapidly decreases [35,36] and ( is constant on large scales, 

i 

with a strongly scale-dependent spectrum oc k~ 2 . The Bardeen potential related to 

i 

( by Eq. (65), grows rapidly and has a flat spectrum [28,8,37] VI oc k°. But equation (55) 
shows that it is only the comoving shear that affects (, and the shear is related to spatial 
gradients of the Bardeen potential, Eq. (63). A scale-invariant Bardeen potential ($ oc k°) 
corresponds to a strongly tilted blue spectrum for the shear (<r oc k 2 ). 

2. Non-adiabatic perturbations 

If we relax the requirement that the matter perturbations are adiabatic (but still assum- 
ing no anisotropic stress) then we no longer have a closed system of equations for ( and $. 
However, we can still estimate $ in the long-wavelength regime given ( and integrating 

2 -H-^+ 5 -±^$^-(l + w)C. (70) 

We also need an initial condition a few Hubble times after horizon exit during slow-roll 
inflation. Starting with the vacuum fluctuation, direct calculation shows that ( at this stage 
is either practically constant (single-field inflation) or only varying slowly on the Hubble 
timescale (multi-field inflation). Through Eq. (70) this gives $ ~ — (3/2) (1 + w)(, and 
hence |$| < |(|. 

We are now going to argue that Eq. (70) will keep |$| < \(\ throughout the super- 
horizon era. A rough argument is the following. Suppose that instead |$| ^> \(\ in some 
super-horizon regime. Then Eq. (70) becomes 

V 1 (ln$)--^f^<-(2/3), (71) 

where we used the energy condition w > — 1 which is always satisfied in scalar field theory 
with a positive kinetic energy. This equation shows that |$| would always be decreasing in 
any expanding universe where |<&| ^> \(\, suggesting that such a regime cannot actually be 
reached starting from the initial condition |$| \(\. 

A more direct argument is to integrate Eq. (70), giving 

F$ = -§ f laa (l + w(a'))F(a')aa')d(lna'), (72) 

where 



InF = 




Assuming that ( is never very much bigger than its primordial value, this will give |<3>| < |C| 
for any reasonable behaviour of w(a). 
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VI. CONCLUSIONS 



In this paper we have shown how local conservation laws (e.g., energy conservation or 
baryon number conservation) can lead to conserved perturbations in cosmology. Whenever 
we have a local continuity equation of the form given in equation (53), then we can construct 
a cosmological perturbation which is conserved after uniform expansion along comoving 
worldlines. 

In particular we have shown that in linear perturbation theory the integrated expansion 
along comoving worldlines between spatially flat slices is just given by the comoving shear. 
Thus on sufficiently large scales (where the shear is negligible) the quantity Xf defined in 
Eq. (54) derived from the conservation equation (53) is conserved. The choice of spatially 
flat slices gives a gauge-invariant definition of the conserved quantity. This is a purely 
geometrical result, whose derivation does not require any gravitational field equations. We 
only require the gravitational field equations in order to estimate the actual comoving shear, 
finding it to be negligible on super-horizon scales. 

The best known example is the curvature perturbation ( = X p , which specifies the total 
density perturbation on spatially flat slices or equivalently the curvature perturbation on 
uniform-density slices. ( is constant on sufficiently large scales (where the comoving shear 
is negligible) for adiabatic density perturbations, for which the local pressure is a unique 
function of the local density and hence the total energy conservation is of the form required 
in equation (53). We have shown that on super- horizon scales, ( coincides with the comoving 
curvature perturbation. 

It is also possible to construct other perturbed quantities, such as the separate curvature 
perturbation Q = X p . for any perfect fluid whose energy is separately conserved [14], or 
d = X m , for any conserved number density, n i? obeying a local conservation equation of the 
form hi = —3rii [23]. 

We also give an expression for the conserved quantity to second-order in the density 
perturbation which may be employed in calculations of the primordial non-Gaussianity of 
density perturbations produced from inflation. 
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APPENDIX 

The Sasaki-Stewart expression for ( 

On super-horizon scales, Eq. (50) becomes 

56 = -3tp . (74) 
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This is valid in any gauge. Choose now a gauge whose slicing is flat at time ti, and uniform- 
density at time t. Integrating from t\ to t, and using ( = —ip on the slice at t, we find that 
on super-horizon scales 

C(x,t) = i f Odt = 6N ss (x,t) (75) 

O Jti 

where JVss(x, t±,t) is the integrated expansion from the flat slice at time ti, to the uniform- 
density slice at time t. This is the expression of Sasaki and Stewart [38], used by them 
to calculate the curvature perturbation at the end of multi-field inflation. It is practically 
independent of the threading that defines the expansion, by virtue of the fact that we are 
dealing with super-horizon scales. 

Our expression Eq. (9) reads ( = —5N. In contrast with the Sasaki-Stewart expression, 
this one is valid only during an era when ( is constant, corresponding to an adiabatic 
pressure perturbation. To understand the relation with the Sasaki-Stewart expression, we 
can integrate Eq. (5) from a uniform-density slice at time t to a flat slice at time ti, both times 
being within the era when the pressure perturbation is adiabatic. Using 3( = Sp/(p + P) 
on the slice at t\ we get the time-independent result 

( = -5N, (76) 

where N = — iVss is the integrated expansion from t to t±. We see that the 'integration 
constant' SN introduced in Section II can be interpreted as a perturbation in the integrated 
expansion between two slices, and that it is equal (as it must be) to — 5N$s- 



Uniform Hubble parameter on comoving slices 

In a given gauge, the perturbation in the expansion 6 with respect to coordinate time is 
given by Eq. (50), which for the flat slicing becomes Eq. (51), 

56 = -a. (77) 

On super-horizon scales this gives \50\/H <C 1, valid for any threading. In other words, the 
expansion with respect to coordinate time is practically unperturbed on flat slices. 

We could instead consider the perturbation in the expansion with respect to proper time, 
given by Eq. (47). On the the comoving slicing one has in the absence of anisotropic stress 
(Eqs. (5.19) to (5.21) of Bardeen [15]; see also [39,33]) 

1Z = -^ com = HA com -H^pj , (78) 

and therefore 

(S6) com = -a. (79) 

Like Eq. (77), this expression is valid for any choice of the threading that defines the ex- 
pansion, and on super- horizon scales it gives \S9\/H <^ 1 practically independently of the 
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threading. In other words, the expansion with respect to proper time is practically un- 
perturbed on comoving slices, for any choice of the threading. In particular the comoving 
expansion is practically unperturbed on such slices, 5H/H <C 1. 

We argued in Section V that <3> is finite on super- horizon scales (in fact, that it is most 
of order the curvature perturbation), and from Eq. (62) this implies that on comoving slices 
the density contrast on super-horizon scales is practically zero, 

\Spcom/p\ < 1 • (80) 

To summarise: on comoving slices, the perturbations in the locally-defined Hubble pa- 
rameter and in the energy density are both negligible in the super-horizon regime. The only 
significant perturbations on comoving slices are therefore curvature perturbation TZ, and the 
pressure perturbation if it is not adiabatic. The statements of the previous paragraph re- 
main true if we replace the comoving slicing by the uniform- density slicing, since we argued 
in Section V that these slicings practically coincide on super-horizon scales. 



The adiabatic condition on the pressure perturbation 

In the text we defined the adiabatic condition on the pressure perturbation as the con- 
dition that the local pressure is a practically unique function of the local energy density. 
Taking it to be absolutely unique, we obtain the familiar adiabatic condition 

8P=(P/p)6p. (81) 

However, on the comoving slicing where Sp is anomalously small, and on the uniform- density 
slicing where it vanishes, it is too strong to require that this expression is valid; there is no 
reason why the slices of uniform pressure should exactly coincide with the slices of uniform 
energy density even if the local pressure is a 'practically' unique function of the local energy 
density. 

An example is provided by single-field slow-roll inflation. During slow-roll the locally- 
defined inflaton field is a practically unique function of proper time, 4>(j), up to the choice of 
origin for r. On super-horizon scales, where spatial gradients are practically negligible, this 
gives practically unique functions p(r) and P(t), making P a practically unique function 
of p. In other words, the adiabatic condition for the pressure perturbation is satisfied on 
super-horizon scales during single-field slow-roll inflation (and even afterwards provided that 
no other field plays a significant role). However, on comoving slices the potential V((f>) is 
uniform, and as a result 

3-Pcom = <5Pcom • (^2) 

This is not in accordance with the strict definition Eq. (81) of an isocurvature pressure 
perturbation. In particular, during slow-roll inflaton p ~ — P(~ V) which means that the 
adiabatic condition for the pressure perturbation is 

5P ~ -5p . (83) 

On a generic slicing this is well-satisfied, but on the comoving slicing it is at variance with 
Eq. (82). All that matters, though, is that both the pressure perturbation and the energy 
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density perturbation are both very small on the comoving slices. Or, to put it differently, 
that the pressure perturbation is very small on uniform-density slices. This is enough to 
ensure that the local pressure is a practically unique function of the energy density, leading 
to the conclusion that ( is constant during single-field inflation. 
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